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SUMMARY 

An analytical method for including a variable 
i sent ro pic exponent in the method of characteristics 
as applied to the design of rocket nozzles has been 
obtained , which is independent of the manner in 
which the initial boundary conditions are specified 
for the problem . The method incorporates existing 
thermodynamic data for either frozen or equilibrium 
composition in the characteristic flow equations of 
s upersonic flo w. 

The design of several nozzles was attempted by 
specification of the initial boundary conditions along 
the nozzle centerline . It was extremely difficult, 
however , to specify the flow distribution along the 
nozzle axis and to obtain solutions to the characteris- 
tic equations yielding short-length nozzles. Selection 
of adequate initial boundary conditions appears to be 
a trial -and -error process based upon previous 
experience . For the case where nozzle length is not 
important, designs are readily obtainable by utilizing 
a geometric scaling procedure . Several bell-shaped 
nozzles have been designed by specifying the initial 
boundary conditions along the nozzle centerline and 
by utilizing geometric scaling. 

A comparison is presented of nozzle contours 
obtained by assuming either frozen or equilibrium 
composition. /I comparison of contours for both 
compositions with either constant or variable isen- 
tropic exponent is also presented. A mmonia, -oxygen 
was used as the propellant combination in the com- 
parison. Significant differences among the con- 
tours of the examples indicated that the assumption 
of constant isentropic exponent is inadequate and 
that actual thermodynamic data should be used in 
the solution of the characteristic flow equations, 

A computation of the variation of vacuum specific 
impulse with axial length showed that considerable 
nozzle length , and hence weight , can be eliminated 
without serious thrust penalties in nozzle designs 


that gradually expand the flow to uniform err it 
conditions. For the example computed \ a reduction 
in axial length of 50 percent resulted in only a 1.6- 
percent reduction in vacuum specific impulse. 

INTRODUCTION 

Generally, the technique described for shock-free 
nozzle design is the method of characteristics, 
which is a linearized solution of the supersonic 
flow equations. Previous application of the 
method of characteristics to nozzle design has 
considered the isentropic exponent 7 to be constant 
throughout the expansion process. (All symbols 
are defined in appendix A.) Reference 1 presents 
a rocket-nozzle design using the method of charac- 
teristics where the isentropic exponent is held 
constant. 

The nature of the expansion process in a nozzle 
is dependent on the chemical kinetics of the gas 
species. Unfortunately, little information exists 
on the chemical kinetics as applied to a rocket 
nozzle, and hence an exact determination of the 
expansion process is impossible. Thermodynamic 
data do exist, however, for the two limiting expan- 
sion processes where the reaction rates are con- 
sidered to be either infinitely slow or infinitely fast. 
For these two processes, the chemical composition 
remains either constant (frozen) or in equilibrium 
throughout the expansion. 

A comparison of one-dimensional isen tropic- 
expansion computations for the frozen and equilib- 
rium cases shows that there is a considerable 
(lifTerence in the area-ratio requirement for a given 
expansion pressure ratio. Also, for both equilib- 
rium and frozen expansion, 7 does no! remain 
constant throughout the expansion process. Con- 
sequently, in computing the contours of nozzles for 
higli-temperature gases, the method of charac- 
teristics should include 7 as a variable. 
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This report contains a description of the analyti- 
cal method for including 7 as a variable in the 
method of characteristics. The analytical method 
uses existing thermodynamic data for frozen or 
equilibrium composition and incorporates these 
data in the characteristic equations of supersonic 
flow. The method for including variable 7 is not 
dependent on the manner in which the initial 
boundary conditions are specified for the problem. 

Included in this report are several bell-shaped 
nozzles that have been designed by the method 
of characteristics for use with the stoichiometric 
combustion products of ammonia and oxygen. 
These nozzles were designed for use in a rocket- 
powered supersonic tunnel described in reference 2. 

A comparison is presented of the nozzle contours 
obtained by using frozen or equilibrium composi- 
tion. A comparison of contours is also presented 
for both compositions with both constant and 
variable isent topic exponent. Also included in 
the nozzle-design technique is a computational 
method for determining the specific impulse at 
each axial station along the nozzle axis. This type 
of computation is particularly applicable to studies 
of rocket nozzles where information on the relation 
between thrust and length is important. From 
the nozzle designs presented herein, some idea 
of the thrust penalties incurred because of shorten- 
ing the length of a characteristic nozzle can he 
obtained. 

METHOD OF ANALYSIS 

AXISYMMETRIC POTENTIAL-FLOW CHARACTERISTIC 
EQUATIONS 

The axisymmetric characteristic equations for 
stead \~, irrotalional, isentropic flow used in this 
analysis were those presented by A. Fcrri in 
reference 3. The characteristic equations are 
obtained by combining the equations of motion 
and continuity. Another relation is derived from 
the geometric properties of the characteristic line. 
The resulting equations for the two characteristic 
lines associated with each point in the How field 
arc 


and, second family: 
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The Mach angle m in equations (1) to ( 4 ) is de- 
pendent on 7 and is therefore a variable at every 
point in the flow field. In figure 1 the velocity 



Figure t. Characteristic lines and velocity vectors for 
a general point in the supersonic stream. 


component and angle relations arc given for a 
characteristic point in the supersonic flow field. 
The same notation used by A. Ferri in reference 
3 has been followed. 


First family: 


FINITE-DIFFERENCE SOLUTION OF CHARACTERISTIC 
EQUATIONS 




0 ) 


C-ttn » i'-0 (2) 

V cos ( 0 +m) V 


Equations (1) to ( 4 ) are numerically integrated 
by rewriting the differential equations in the form 
of finite differences and following a point-to-point 
calculation of the flow properties throughout the 
flow field. Figure 2 shows two points, A and fi, 
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Figure 2. — Point calculation system in the characteristic 
network. 


the differential equations (2) and (4), so that they 
become 
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in the characteristic network, whoso properties 
are known. If the distance between character- 
istic points is kept small, the tangents of the 
characteristics can be assumed to represent the 
characteristics themselves for a first approxima- 
tion. The tangents of the characteristics from 
two different families at points 7? and A, that is, 
tan (Sff—fjLn) and tan (0 A + ju A ), respectively, inter- 
sect at point F. The flow properties and location 
at Fare determined by first solving equations (1) 
and (3) simultaneously for x c and ?/ c ; and second, 
knowing x c and y c , equations (2) and (4) are solved 
simultaneously for V c and B c . Thus for the first 
approximation, the location of point Fis given by 
equations (1) and (3) in finite-difference form: 


By solving for V c by elimination of 6 C} the follow- 
ing relation is obtained: 
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and solving for the 

flow angle d Cj tl 

ie following 

relation is obtained: 
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Solving for x c from equations (5) and (6) yields 

_ JjA—ya+*B tan (B b —pb)—Xa f an ( B a +Pa) ( 7 

Xc tan (B b — ix B ) tail (B a +Va) 

and for y Ct 

Vc = y . 4 + (x c — x A ) tail (6 A + ix A ) (8) 

Likewise, finite-difference equations can replace 


The values obtained from equations (7), (8), (II), 
and (12) are only a first approximation, since all 
the coefficients in the equations were based on the 
properties at points A and B. Usually a second 
iteration of the equation is required where the 
coefficients of the differential equations are aver- 
aged between .1 and F, and B and F. Thus, the 
second iteration point is located at F ; , and equa- 
tions (7) and (8) become 


y A — 2/w+2 t* an gfl)+tan (B c — ix c )]x B ~ [tan (B A +jx A )+iim (B c +^c)]7 A 

^ [tan (B b — ///,)+ tan (B c — ju c ) — tan (0 a + Ma)“ tan(0 c +/*c)] 
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and 
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Similarly, equations (11) and (12) become 
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Apparently, in order to iterate the equations 
after the first approximation, the value of Much 
angle ji at C must be obtained so that the coeffi- 
cients can be averaged. Tn order to calculate 
the Mach angle at C } the speed of sound must 
be determined from the thermodynamic proper- 
ties of the gas. If the properties change signifi- 
cantly between the first and second iterations, a 
third iteration may be required. After the flow 
properties at have been determined within the 
desired accuracy, point C 9 becomes a starting 
point for the calculation of future points of the 
network. 


THERMODYNAMIC ANALYSIS 

The isen tropic exponent is an important param- 
eter in all compressible-flow processes but can be 
considered to bo constant in any low-temperature 
airflow analysis. However, in a rocket- nozzle 
operating with high-pressure and high-temperature 
gases, the dissociation effects of the gas should 
be included in the supersonic characteristic 
flow solutions. The dissociation and reasso- 
ciation phenomena and chemical kinetics cause 
the molecular weight of the combustion gases to 
vary with both time and position in the nozzle. 
When the molecular weight varies, the ison- 
tropie exponent (which enters into the computa- 
tion of the speed of sound) cannot be expressed as 
the ratio of specific beats, but must include a 
function of the rate of change of molecular weight 
with pressure. The expression for the insen- 
tropic exponent can easily he derived from its 


definition by using the tables of reference 4 and 
the equation of state. Tlenco, 



03 ) 


The relation between molecular weight and 
pressure for evaluation of the isen tropic exponent 
is dependent on the chemical kinetics of the 
chemical constituents of the reacting combustion 
gas. Unfortunately, the necessary chemical kinet- 
ics arc not known for current rocket propellants 
except for the two limiting cases of frozen and 
equilibrium composition. A generalized method 
for computing the thermodynamic properties of 
a combustion gas assuming equilibrium or frozen 
composition is given in reference 5, and data for 
a number of rocked propellant combinations are 
presented in references G to 10. 

The dal a for frozen or equilibrium composition 
can be incorporated into the finite-difference 
solution of the characteristic equations by utilizing 
the appropriate thermodynamic relations. In- 
serting the thermodynamic data into the charac- 
teristics solution can be achieved with only minor 
additions to the method of characteristics. When 
high-speed automatic-computing machines are 
utilized to obtain a characteristic solution, adding 
the thermodynamic data to the solution of the 
method of characteristics does not add undue 
complications. 

Any thermodynamic property can always he 
considered to be only a function of any two other 
thermodynamic properties. For an isen tropic 
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process it is therefore possible to express all 
thermodynamic data for a given entropy level in 
terms of only one thermodynamic property. If 
the enthalpy is considered to be some function 
<p of the natural logarithm of the pressure, then 
from the energy equation for isen tropic flow the 
static-temperature molecular-weight ratio can 
be obtained from 


r T_ (UI _ d<p 

in d In p din/? 


(14) 


Similarly, the following y function can be obtained 
by using the tables of reference 4 and the equation 
of state: 


T 

7—1 ^ dV j d(p 

7 d In p d(ln p) 2 j d(ln p) 


(15) 


Equation (15) gives the iscntropic exponent 
based on the thermodynamic properties of the 
rocket gases for each point in the nozzle. If the 
local-temperature -molecular-weight ratio and the 
local iscntropic exponent are known, the speed of 
sound at the point in the flow is determined from 


a= 



(16) 


The thermodynamic properties are related to 
the velocity by the energy equation in the form 


V— K q (TI c — 77) 1 72 — K q (77 c — <p) 1 12 (17) 


Equations (14) to (17) give the relations neces- 
sary to evaluate the Mach angle n for any given 
velocity. 

The enthalpy function <p was obtained by 
fitting the data for frozen or equilibrium com- 
position over the pressure range for the nozzle 
design. The following numerical relation using 
In p as a parameter was fitted to the thermo- 
dynamic data: 


<P = TT= Tv, + Ki (1 1 1 p)+K,(\n p) 2 

+K&n p) 3 +Iu ( In pY+K&n p)* (18) 


The preceding method for including variable 7 in 
the method of characteristics is general and is 
independent of the manner in which the initial 
boundary conditions are specified for the problem. 


APPLICATION OF THE CHARACTERISTIC EQUATIONS 

The method for incorporating the thermo- 
dynamic data into the characteristic equations 
was used to design a shock-free divergent section 
for the ammonia-oxygen rocket tunnel with a 
chamber pressure of COO pounds per square inch 
reported in reference 2 . For design of nozzles 
for tunnel application it is necessary to specify 
the exit flow and then determine the contour that 
will produce that flow. A weight flow of 19.06 
pounds per second and an exit pressure of approxi- 
mately 1 .7 pounds per square inch absolute were 
desired for operation at stoichiometric conditions. 
This corresponds to an exit Mach number of 
approximately 4.0. The thermodynamic data 
incorporated in the characteristic calculations are 
reported in reference 10 . The pressure distri- 
bution along the nozzle centerline from throat to 
nozzle exit was prescribed, and the corresponding 
properties of the gas were determined from 
equations (14) to (18). (The flow angle along 
the centerline of an axisymmetric nozzle is zero.) 
Thus, the properties at any two adjacent points 
along the centerline of the characteristic network, 
such as D and E of figure 2 , provide the initial 
information for the solution of the characteristic 
equations (5), ( 6 ), ( 11 ), and ( 12 ). The finite- 
difTercncc computation proceeds outward from 
the nozzle centerline. The details of the com- 
putational procedure are given in appendix B. 

Solutions to the characteristic equations were 
obtained for centerline pressure distributions 
defined as cosine, parabolic, and cubic functions 
of the axial length. The functions, shown in 
figure 3, were defined by 

1 7T.T 

Cosine: In p— ^ (In p ,— In p e ) cos-^- 

+1 (In /b+ln p e ) 
(19) 

Parabolic: In p=l~ 2 (x— Z) 2 (ln p t — In p e ) 

+lr ij>. (20) 

C„1A: In 

( 21 ) 
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Ficthre 3. — Cen tori inn pressure input distributions. 

The expansion length l chosen was 12 inches, 
and the exit pressure was 1.7 pounds per square 
inch absolute. The pressure at the throat was 
chosen as the value obtained by using equilibrium 
composition. The coefficients of equations (19) 
to (21) were determined from the following as- 
sumed boundary conditions: 


at x—x t ~ 0, 

d In p/dx= q, 

P=P 

at z=l— 12, 

d In p/dx=0, 

p—p< 


In addition to varying the centerline pressure 
distribution using equilibrium composition, a 
comparison was made of the effect of frozen and 
equilibrium composition upon the solution with 
a fixed centerline pressure distribution. A com- 
parison was also made (for each composition) of 
the effects of both fixed and variable isen tropic 
exponent upon the solution when using a fixed 
centerline pressure distribution. 

THROAT CONDITIONS 

An initial value problem is defined when the 
centerline pressure distribution for the supersonic 
portion of the nozzle is specified. The solution of 
this initial value problem does not encompass the 
entire supersonic portion of the nozzle, as shown 
in figure 4. No solution is obtained for the region 
between the sonic line arid the characteristic 
that passes through the sonic point on the center- 
line. For the case where the throat curvature is 
very small, the undefined region of the supersonic 
flow field becomes small, and the flow in this region 



Ficube 4. Sketch of undefined region of the character- 
istic solution with boundary conditions specified along 
nozzle centerline. 


may be neglected. This is the case for the nozzle 
designed for the rocket tunnel of reference 2. 

ADDITIONAL CALCULATIONS 

Weight flow. — Any streamline determined from 
the results of the characteristic solution generates 
a shock-free rocket-nozzle contour. A weight- 
flow computation was carried along with the 
characteristic calculations so that the data arc 
available for streamline calculations. For in- 
stance, consider the weight flow across the char- 
acteristic line AC of figure 2: 

w c —w A =j ) (pu2irydy—2wpvy dx) 
which reduces to 

w c -w A = 2r J V[cos g- ydy (22) 

A numerical integration of equation (22) be- 
tween characteristic points gives a weight-flow 
history of the characteristic network, and the 
streamlines arc easily obtained by interpolation. 

Vacuum specific impulse. The vacuum- 
specific-impulse (assuming the flow discharges to a 
vacuum) calculation can bo executed as the char- 
acteristic flow network is being developed. This 
computation is of considerable importance when 
the contours of rocket nozzles are being computed 
because it indicates the effect of nozzle size on 
thrust. For these calculations the following equa- 
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lion was used along the family of AC character- 
istics: 

V A cos 0 A -\-V c cos 6 C \ 

""“ 2 " / 

(23) 

where all values of l, w, and F arc integrated 
values from the centerline to the point repre- 
sented by the subscript. 

RESULTS AND DISCUSSION 

The characteristic calculation using equilibrium- 
composition data with the cosine-function pressure 
distribution as input showed a coalescing of 
one family of characteristic lines a short distance 
downstream of the throat at a weight flow of 
about 2 pounds per second. This coalescing indi- 
cates that it is impossible to obtain a solution 
having the desired flow rate for this particular 
centerline pressure distribution. 

To learn the effects on the characteristic net- 
work of changing the centerline pressure distribu- 
tion, the parabolic and cubic inputs of figure 3 
were also tried. The results were essentially the 
same as experienced with the cosine input; namely, 
a coalescing of one of the families of characteristic 
lines at a weight flow of approximately 2 pounds 
per second. Thus, it appears that, if the flow- 
expansion length is too short, it will be extremely 
difficult to find a characteristic solution with a 
pressure distribution along the centerline being 
described by any function of the family of the 
simple functional relations considered herein. 
'Reference 1 1 gives a somewhat detailed discussion 
of t lie selection of centerline pressure distributions 
for the related two-dimensional case and points 
out this problem of coalescing of characteristics. 
Tt is indicated that only a trial-and-error method 
based on previous experience can give the desired 
characteristic solution. Although no solution 
was attainable with the inputs of equations (19) 
to (21), it was possible to scale the solution using 
the parabolic input to obtain a shock-free contour 
for the rocket-tunnel weight-flow requirements. 

SCALING 

Examination of equations (1) to (4) shows that 
all the coordinates / and y appear only as ratios. 
Tt is therefore feasible to use a geometric scaling 

527411 GO 2 




procedure to obtain nozzle contours with weight 
flows different from those of existing solutions of 
the characteristic equations. Since the radial 
dimension of the nozzle is proportional to the 
square root of the weight flow, the axial dimension 
must also vary as the square root of the weight 
flow in order to preserve geometric similarity. 
For the streamline with a weight flow of 2.1 pounds 
per second obtained with the parabolic centerline 
pressure distribution, a sealing factor of 3 was 
necessary to scale to a nozzle contour with a weight 
flow of 19 pounds per second. 

The centerline pressure distribution is alt (‘red 
when geometric scaling is used to obtain a new 
solution. When a 2.1 -pound-pcr-sccond weight- 
flow streamline is scaled up to the nozzle with a 
1 9-pound-per-sccond weight-flow rate, the pres- 
sure gradient along the centerline is reduced by a 
factor of 3. The sealed parabolic centerline pres- 
sure distribution is shown in figure 5. 



0 4 8 12 16 20 24 28 32 36 


Axia! nozzle length t x } in. 

Ficvre 5. Parabolic pressure input scaled three times. 

Nozzle designs based upon equilibrium and 
frozen composition with both variable and con- 
stant isentropic exponent for each composition. 
The thermodynamic properties of the high- 
temperature combustion gases in a nozzle are usu- 
ally somewhere between the extremes determined 
by assuming equilibrium or frozen composition. 
The effect of real thermodynamic data on nozzle 
design can be indicated, however, by examining 
the effects of the assumption of either frozen or 
equilibrium composition. The effects of con- 
sidering a constant y can be examined by compari- 
son with solutions with variable y. 

The streamlines resulting from the character- 
istic calculations using the input of figure 5 are 
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shown in figure G. The flows obtained by assum- 
ing a variable 7 and either frozen or equilibrium 
composition are presented. Also shown in figure 
G are (he isobars and streamlines for both frozen 
and equilibrium composition when the ison tropic 
exponent is held constant. The isentropic ex- 
ponents used in the two constant 7 solutions were 
the values corresponding to each chamber condi- 
tion. A value of 1.1402 was used for (lit* isen- 
tropic exponent when equilibrium composition 
was used, while a value of 1.2054 was used with 
the frozen composition. All other thermody- 
namic properties used in the two constant 7 solu- 
tions were identical with those used in the cor- 
responding variable 7 solutions. 

As shown in the four examples of figure G, the 
major portions of the axial length are spent in 
turning the flow very slowly to the axial direction. 
This is typical of nozzles designed for tunnel appli- 
cation where gradual accelerations and uniform 
exit flow are desired. 

A comparison of the 4- and 18-pound-pcr-second 
streamlines for frozen and equilibrium composition 
with both fixed and variable 7 for each composi- 
tion is given in figure 7. There is a significant 
difference in the flow contours for the two vari- 
able 7 solutions. For 1 the examples computed, the 
solution with equilibrium composition gives an 
area in the axial-flow section approximately 20 
percent greater than the solution with frozen com- 
position. For frozen composition, the streamlines 
obt ained by assuming a constant 7 closely approxi- 
mate the streamlines obtained by using a variable 
7 . A slightly larger area ratio results from the 
constant-isentropic-exponcnt solution in the mid- 
dle portion of the nozzle; however, the exit areas 
are nearly equivalent. 

The solution based on constant isentropic ex- 
ponent with equilibrium composition is significantly 
different from the oilier three solutions. The 
characteristics of the solution with equilibrium 
composition and a constant 7 coalesce at a weight 
flow of approximately 12 pounds per second, and 
hence no solution is possible for greater weight 
flows. Appreciable differences in the streamlines 
at lower flows are noted. At a weight flow of 4 
pounds per second, the differences in streamline 
shape in the middle portion of the contour between 
the two equilibrium-composition cases are as great 
as the differences between the two variable 7 cases. 

The differences among the solutions based on the 


different assumptions show that the variation of 7 
is an important factor in the design of nozzles oper- 
ating at high temperatures. For more-reactive 
propellant combinations, such as hydrogen-fluorine, 
there is an increase in the difference in area-ratio 
requirements for frozen and equilibrium composi- 
tion. In addition, the variation of 7 between 
throat conditions arid nozzle exit conditions in- 
creases for both frozen and equilibrium compo- 
sition. These effects can be seen by comparing 
the onc-dirncnsionul isentropic data of references 
G to 10 . For more-reactive propellant combina- 
tions it would therefore become more important 
to include thermodynamic data that comprehend 
the effect of a variable isentropic exponent. 

Vacuum specific impulse.— The variation of 
vacuum specific impulse with axial length is of 
great interest in the design of rocket nozzles. 
The results of the computation of vacuum specific 
impulse for the solution obtained with constant 

7 and equilibrium composition are given in figure 

8 for the 4-pound-per-second streamline. Con- 
siderable nozzle length, and hence weight, can be 
eliminated without serious thrust penalties. If 
the axial length of the 4-pound-per-sccond stream- 
line is reduced by 50 percent, there is only a 1.6- 
percent loss in vacuum specific impulse. As the 
length is further reduced, the loss in vacuum 
specific impulse increases rapidly. A reduction in 
length of 75 percent will give a 9.7-percent reduc- 
tion in vacuum specific impulse. It is hence 
practical to eliminate a large portion of the nozzle 
that has the greatest surface area per unit axial 
length and thus greatly reduce the weight without 
serious penalties in vacuum specific impulse. 

SUMMARY OF RESULTS 

An analytical method for including a variable 7 
in the method of characteristics as applied to the 
design of rocket nozzles lias been obtained. The 
method of including a variable 7 is independent of 
the manner in which the initial boundary condi- 
tions are specified for the problem. The method 
incorporates existing thermodynamic data for 
either frozen or equilibrium composition in the 
characteristic flow equations of supersonic flow. 

The design of several nozzles was attempted by 
specification of the initial boundary conditions 
along the nozzle centerline. It was extremely 
difficult, however, to specify the flow distribution 
along the nozzle axis and to obtain solutions to the 






i cure 6. — Continued. Streamlines and pressure contours for nozzles with a parabolic centerline pressure distribution. 
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FiorRE 7. — Effect of frozen and equilibrium composition and isentropic exponent on nozzle contour for parabolic 
cenlerline pressure distribution and stoichiometric mixtures of ammonia and oxj’gcn. 
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Axial length, x t in. 


Figure 8* — Variation of vacuum specific impulse with 
axial length for a streamline. Weight flow, 4 pounds 
per second; equilibrium composition; constant isentropic 
exponent. 

elm met eristic equations yielding short -length noz- 
zles. Selection of adequate initial boundary 
conditions appears to be a trial-and-error process 
based upon previous experience. For the case 
where nozzle length is not important, designs are 
readily obtainable by utilizing a geometric scaling 
procedure. Several hell-shaped nozzles have been 
designed by specifying the initial boundary condi- 


tions along the nozzle centerline and by utilizing 
geometric sealing, 

A comparison of nozzle contours obtained by 
assuming either frozen or equilibrium composition 
is presented. A comparison of contours for both 
compositions with either constant or variable 
isentropic exponent is also presented. Ammonia- 
oxygen was used as the propellant combination in 
the comparison. Significant differences among the 
contours of the examples indicated that the as- 
sumption of constant y is inadequate and that 
actual thermodynamic data should be used in the 
solution of tlu* characteristic flow equations. 

A computation of the variation of vacuum 
specific impulse with axial length showed that 
considerable nozzle length, and hence weight, can 
be eliminated without serious thrust penalties 
in nozzle designs that- gradually expand the flow 
to uniform exit conditions. For the example 
computed, a reduction in axial length of 50 percent 
resulted in only a l.G-perccnt reduction in vacuum 
specific impulse. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February It, 1050 


APPENDIX A 


SYMBOLS 


A,B,C,D,E 

a 

c v 

V 

j 

g 

II 

I 


, K q 

^ 5j^6 

L 


l 


M 

.// 

ra 

P 

4 


characteristic point identification 
local speed of sound 
specific heat at constant pressure 
specific heat at constant volume 
integrated thrust from centerline 
to point represented by subscript 
arbitrary factor, eqs. (Bl) and 
(B2) 

acceleration of gravity 
enthalpy 

integrated specific impulse in vac- 
uum from centerline to point 
represented by subscript 
constant for conversion of units 
thermodynamic data curve-fit ting 
coefficients, eq. (18) 
characteristic line 
coefficient in characteristic equa- 
tion 

flow expansion length, axial dis- 
tance along centerline for flow to 
expand from p t to p e , /— (length 
of supersonic portion of nozzle) 
— yjt an p e 
Mach number 

coefficient in characteristic equa- 
tion 

molecular weight 
local static pressure 
constant in eq. (21) 


R 

T 

t 

u 


V 


v 

10 


X 

V 

y 

o 

p 

p 

<p 

Subscripts: 

a,b,c,d 

C 

e 

cst 

3 

T 

t 

Superscript: 


universal gas constant 
local static temperature 
tangent to characteristic line 
component of velocity in x-dircc- 
lion 

flow velocity 

component of velocity in ^/-direc- 
tion 

integrated weight flow from center- 
line to point represented by 
subscript 

axial distance along nozzle 
radial dimension of nozzle 
isentropic exponent, 

7 = (d In pjb In p) s 
flow angle 

Mach angle, g — sin 1 1/71/ 
density 

enthalpy function 

characteristic point identification 

chamber 

exit 

estimated 
constant entropy 
constant temperature 
throat 


second iteration 



APPENDIX B 

CALCULATION PROCEDURE 


Unless extremely coarse grids are used, the 
number of equations and iterations involved in a 
nozzle characteristic network compulation make 
hand-eompulat ional methods formidable. The 
method of characteristics applied to nozzle design 
as presented in this report therefore has been pro- 
grammed and computed on an IBM 053 digital 
computer. Essentially the following procedure 
was employed in programming the method: 

The calculation using the method of charac- 
teristics is initiated in the supersonic portion of 
the nozzle using an assumed pressure distribution 
along the centerline as the initial boundary con- 
ditions. The centerline static -pressure distribu- 
tion is specified from the nozzle throat to the exit 
in terms of the natural logarithm of the static 
pressure as a function of the axial length along 
the nozzle centerline. The method of character- 
istics imposes no restriction on the centerline 
pressure distribution except that the function and 
at least the first derivative of the function should 
be continuous. Specification of a centerline pres- 
sure distribution will uniquely determine a flow 
field if a solution exists. 

After the centerline pressure distribution is 
specified, grid points for the characteristic net- 
work are selected along the nozzle axis. A max- 
imum grid spacing should be used; however, the 
spacing should also be fine enough to insure the 
accuracy desired in the finite-difference solution. 
A few trial grid spacings will determine the proper 
interval. 

For convenience, the computation may be ini- 
tiated at either the exit or the throat of the noz- 


zle. The properties of the fluid at each of two 
adjacent points, D and E, along the nozzle axis 
may be determined from equations (14) to (18). 
Having a complete description of the flow along 
the nozzle axis, the characteristic equations can 
be solved to obtain the properties of the flow at 
the point C. As shown in figure* 2, point C is de- 
fined by the intersection of the 0-f ju characteristics 
originating from point A and the 6 — g character- 
istic originating from point B. 

The d ('termination of the flow properties at 
point C is an iterative process. The pressure and 
flow angle at C are initially estimated by any ar- 
bitrary scheme. The following relations wore 
used to obtain the initial estimates: 

In p c =J(\n Pa — hi p B )+ In p B (Bl) 

and 

e c =(o A -e B )/+e B (B2) 

whore / is any arbitrary factor. 

Based on the estimate of pressure and flow an- 
gle, a complete description of the flow properties 
at point C can be obtained from equations (14) 
to (18), inclusive. 

The characteristic equations (7), (8), (11), and 
(12) can now bo rewritten to provide first the 
location of the point C and then a check on the 
assumed values of flow angle and pressure. All 
coefficients in equations (7), (8), (11), and (12) 
are evaluated using properties that are the aver- 
age of the estimated values at point C and the 
values at either point A or B. Thus, equations 
(7), (8), (U), and (12) become 


= Vb+Xh [ 

r 


tan (6 b — M«)+tan (9 c ,esi— Me, ««) 


tan r 6 A 4- /i, ) -j 1 an (+ , est ~1~ Me, f st 1 
2 


] 

rx 


tan (0B— Ms) + tan (0 (7 ,est ^C,est) 


tan (0,,+^) + tan (0, 


est + Mc,f.sc) r m 


Vc—y. 


, x Ran 

a + ( T c — ?a) I 


( 0 .- 1 +a.i) + tan (0 C , esi + MC, 


est) 


(B4) 
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’4[(F 4 -|-Fc,„i)'' l (1iin /iA + tau He , ««) -1 + (Vn+Vc , „,) ^ 1 ( I )' n Me+tan Mr, «/)"'! 

r 2 r 2V a . ft 1 

l tail MA + fan Hc.en I Fj + Ve en 2 c I 


-Vl>— + ^ d #“ <,c- '- ) ' 

sf At » C,<st ^ 


tan //^+ta 


n tftw — + B ~^~n c,€ - (i^c - ^"V (B5) 

n Me, est L» B\ V C.est * J J 


_ pa c-y 

C, ( St D C.csfT T A 


tan /iA + tan juc.r*/ L^c.««+ 


c— ^4) 


Equation (BO) provides a direct check on tlie 
estimated value of flow angle. Equation (B5) 
provides the velocity at point C, which is an 
indirect check on the assumed static pressure. 
The flow angle from equation (B6) and the velocity 
from equation (B5) must agree with the assumed 
value of flow angle and the value of velocity 
obtained from equation (17) based on the assumed 
value of static pressure. Normally , the values of 
flow angle and velocity will not agree with the 
values initially assumed, and therefore additional 
iterations will be required to obtain convergence. 

An estimate of the necessary change in the 
initial assumptions of flow angle and static pressure 
can be obtained by differentiating equation (17). 
Hence, 

dF -K\ d* . KJ 
d In p 2V d In p Vrn ’ 

The static pressure and flow angle for the second 
and all additional iterations are then given by 

8c' = 0 c +A($ c )/ 2 (B8) 


i npc , =] „ Pc+{T/m ^— A{Vc) 


where AV C and A9 C are the errors in the velocity 
and flow angle, respectively, in the previous 
iteration. This iteration procedure usually con- 
verged in three or four iterations. 

The preceding procedure was carried out for 
each pair of adjacent points along the nozzle 
centerline. Tf there are n points specified along 
the nozzle axis, there will result («-l) points of C f . 
The calculation was then reinitiated using pairs 
of these (n-1) points of C 7 as A and B points, and 


the process was repeated until the entire flow 
field of the nozzle was computed. 
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